Coulomb breakup in a Transformed Harmonic Oscillator Basis 
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The problem of Coulomb breakup in the scattering of a two-body loosely bound projectile by 
a heavy target is addressed. A basis of transformed harmonic oscillator (THO) wave functions is 
used to discretize the projectile continuum and to diagonalize the Hamiltonian of the two-body 
system. Results for the reaction 8 B+ 58 Ni at subcoulomb energies are presented. Comparison 
of different observables with those obtained with the standard Continuum Discretized Coupled- 
Channels (CDCC) method shows good agreement between both approaches. 
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The collision of a weakly bound system with a tar- 
get represents a challenging and interesting problem in 
quantum physics, including atomic, molecular and nu- 
clear physics. A proper understanding of the process re- 
quires an appropriate treatment of the unbound part of 
the spectrum of the loosely bound system. In nuclear 
collisions, the problem was first addressed in the context 
of deuteron scattering. One of the first successful ap- 
proaches to this problem was the Continuum Discretized 
Coupled-Channels (CDCC) method, firstly introduced 
by Rawitscher filLand later developed and employed by 
other groups 011 0]. Within the CDCC method, the 
reaction process of a loosely two-body projectile and a 
heavy structureless target is treated within a three-body 
picture. The idea of the method is to represent the con- 
tinuum part of the two-body projectile spectrum by a 
finite set of square integrable states. To this end, the 
continuum is divided into a finite set of energy inter- 
vals. For each interval, or bin, a representative function 
is constructed by superposition of true scattering wave 
functions within that interval. By construction, the set of 
functions obtained in this way are normalizable and mu- 
tually orthogonal. By projecting the Schrodinger equa- 
tion onto the bound and bin wave functions, a set of 
coupled equations is obtained. The method has been ex- 
tremely successful and is a reliable reference for any other 
alternative method. In spite of this, it has been criti- 
cized by some authors due to the apparent arbitrariness 
in the definition of the continuum bins |5j. Furthermore, 
the generalization of the method to treat the three-body 
continuum is not obvious. 

As an alternative to the binning procedure, other 
methods have been used to represent the continuum spec- 
trum of a two-body system by a discrete and finite set of 
square integrable, or L 2 , states such as Laguerre 0,0 or 
Gaussian H Q functions. A family of these states, usu- 
ally called pseudo states, is used to diagonalize the Hamil- 
tonian of the two-body system. The resulting eigen- 
states are then used within the CDCC calculation in 
exactly the same way as the continuum bins. In this 
context, a new basis suitable for continuum discretiza- 
tion has been proposed and applied to the scattering of 



two-bod y c omposite systems in a series of recent works 
[H El El EH 13 EB The method generates a dis- 
crete representation of the continuum spectrum starting 
from the ground state wave function, which is the only 
needed input. A local scale point transformation (LST) 
[Ifil IrH ] that transforms this function into the harmonic 
oscillator ground state wave function is defined. Once 
the LST is obtained, the inverse transformation produces 
from the harmonic oscillator states the wave functions 
that represent the continuum (and other bound states 
if they exist) of the two-body system. The correspond- 
ing pseudostate set is known as the Transformed Har- 
monic OscillatorTTHO) basis. The method was first de- 
veloped in Ref. [1Q] for simple one dimensional problems 
and later extended in order to check its applicability and 
limitations El El El El El- In particular, in it 
was shown that the combination of the THO discretiza- 
tion method with the coupled channels technique, named 
CDCC-THO, can be useful to describe continuum effects 
in nuclear collisions. In [12], several restrictions were im- 
posed for the sake of simplicity: the deuteron was taken 
to be a pure s-state, the effect of Coulomb breakup was 
neglected, and the nuclear interaction of the proton and 
the neutron with the target was accepted to induce the 
coupling just to s-wave breakup states. Moreover, the 
deuteron ground state was represented by an analytical 
wave function. Within this simplified scenario, it was 
shown that the method was in excellent agreement with 
the standard CDCC method. 

In this work we apply the CDCC-THO method to the 
reaction 8 B+ 58 Ni at 25.8 MeV, in which the 8 B pro- 
jectile is modelled as a weakly bound two-body system, 
proton+ 7 Be. This reaction has been extensively analyzed 
within the standard CDCC formalism 0, E| and the 
CDCC with a Gaussian basis j2jj]. The motivations of 
this work are the following. First we would like to show 
that the method can be equally implemented in the usual 
situation in which the ground state of the system is not 
known analytically but is calculated numerically. Sec- 
ond, we aim to show that the method can be generalized 
in a straightforward way to describe the continuum part 
of the spectrum for arbitrary partial waves. Since the 8 B 
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ground state corresponds mainly to a p— wave configura- 
tion, the THO method provides a discrete representation 
for the 1=1 continuum. We will show below, however, 
that it can be easily extended to generate also a rep- 
resentation for the £ — continuum. Finally, we will 
show that, despite the finite extension and exponential 
asymptotic behaviour of the THO basis, these states are 
suitable to describe situations dominated by long range 
interactions, such as the dipole Coulomb couplings aris- 
ing in collisions of a loosely bound projectile by heavy 
ions. 

For that purpose, we first present a brief review of the 
the basic formulation of the THO method and describe 
how it can be extended to different continuum configura- 
tions. Then, we present calculations, in which the THO 
wave functions are used for a coupled channels calculation 
(CDCC-THO), to show the convergence of the method 
for several magnitudes related to the 8 B+ 58 Ni reaction 
at 25.8 MeV. We use the standard CDCC calculations as 
a reference to compare with. 

The standpoint of the THO method is the ground state 
wave function of the two-body system, here denoted as 
4>b,£ ( r )- The subscript £q represents the inter-cluster 
relative angular momentum. To simplify the notation, 
we do not consider the intrinsic spins of the fragments. 
By use of a local scale transformation (LST) 0, Il7j . 
this wave function is converted into a harmonic oscilla- 
tor (HO) wave function. Thus, the function s(r) defining 
the LST is given by 

Frhi 

Mr) = \j Tr <i>Z?M r ))- (!) 

where <j>^f (s) is the radial part of the HO wave function 
for the orbital angular momentum £q. Once the s(r) func- 
tion has been obtained, the THO basis is generated by 
applying the same LST calculated for the ground state to 
the rest of HO wave functions. Due to the simple analyt- 
ical structure of the harmonic oscillator wave functions, 
this is equivalent to multiply the ground state function 
by the appropriate Laguerre polynomials il +1 ^ 2 (s) 0] 

<f (r) = [ S (r)]^ Ll+ 1 /2 (s(r)2) ^ o(r) . (2) 

Notice that, by construction, the family of functions 
are orthogonal and constitute a complete set. 
Moreover, they decay exponentially at large distances, 
thus ensuring the correct asymptotic behaviour for the 
bound wave functions. However, in general, these func- 
tions are not eigenstates of the internal Hamiltonian. 
Then, the following step is to diagonalize the Hamil- 
tonian using a truncated THO basis. As a result of 
the diagonalization a new set of functions, denoted 
{4> I j,e.{ r )'i j — 0, . . . , N}, with eigenvalues eo, . . . , ejv are 
obtained. Here, N+l is the number of functions retained 
in the THO basis, j = standing for the ground state. 

Thus 0oA,(r) = ^o,^°( r ) = </Vo( r ) and so e = e&, while 
the rest of eigenstates constitute our representation of the 
bound and unbound energy spectrum. Those eigenstates 



at negative energy will represent of bound states of the 
system. In particular, for £ < £ , the diagonalization gen- 
erates negative energy states, which correspond to occu- 
pied states. These states are Pauli forbidden states, that 
should be excluded from the calculations. 

In 0|, the transformation given by Eq. J3J was de- 
fined only for continuum states with the same angular 
momentum I as the ground state. In the present work, 
we generalize this procedure by applying the same trans- 
formation to obtain THO states for t ^ £ . In particular, 
we will apply the method to describe the I = and £ — 1 
continuum states of the 7 Be+p system. 

To simplify the calculations, we neglect the 7 Be and 
valence proton spins. Therefore, the 8 B ground state is 
given by a £ = 1, J = 1 configuration. Using Eq. 
we generate the THO basis for the £ = and £ — 1 
states. Then, we retain a finite number of states and 
diagonalize the 8 B Hamiltonian in this truncated space. 
The resulting eigenvalues and their associated wave func- 
tions constitute a discrete and finite representation of the 
continuum for the 7 Be+p system. 

The bound and continuum states are used to generate 
the diagonal and non-diagonal coupling potentials that 
enter the system of coupled equations. In both the stan- 
dard CDCC and CDCC-THO methods, these coupling 
potentials are defined as 

U a ,a'{R) = (</>n/|£q"p-Ni] + ^[7Be-Ni] 1 4>n' ,£> ) , (3) 

with a = {n,£}, a' — {n',£'}. The internal wave func- 
tions 4> n /(r) are represented by either continuum bins or 
the THO basis, depending on the discretization method. 
The proton- 7 Be nuclear interaction was represented by 
a simple Woods-Saxon form with the potential param- 
eters of Esbensen and Bertsch (U- This interaction is 
used to generate the ground state wave function for the 
8 B system as well as for the diagonalization of the 8 B 
Hamiltonian in the THO basis. The proton- 58 Ni interac- 
tion was taken from the parametrization of Becchetti and 
Greenless [22]. For the 7 Be- 58 Ni system, we used the in- 
teraction of Moroz et al. 2."!| following the choice of other 
authors 0,H3]. All possible couplings (s — s, s — p and 
p — p) were included in the calculations for multipolar- 
ities A=0,1 and 2. In the standard CDCC calculation, 
the continuum spectrum was discretized into energy bins 
of equal momentum width and up to a maximum exci- 
tation energy of e max = 10 MeV. In particular, we used 
N = 16 bins for £ = 1 and N = 32 for I = 0. The 
bin wave functions were calculated up to 100 fm. For a 
meaningful comparison, we take the same potential pa- 
rameters as in the CDCC-THO calculation. The coupled 
equations were solved for projectile-target orbital angu- 
lar momenta up to L max — 500 and integrated up to a 
radius i? ma x = 500 fm. 

In Fig. we represent the eigenvalues obtained upon 
diagonalization of the 8 B Hamiltonian in a THO basis 
with N = 20 states. The vertical axis corresponds to the 
linear momentum ki — y/2jlel/h, where [i is the proton- 
7 Be reduced mass. It should be noticed that these eigen- 
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FIG. 1: Eigenvalues obtained upon diagonalization of the 
internal Hamiltonian in a THO basis with N = 20 states for s 
and p-waves. The vertical axis corresponds to the proton- 7 Be 
linear momentum. The p-state at negative energy corresponds 
to the ground state. 



values tend to concentrate close to zero energy. As the 
excitation energy increases, the energy distribution be- 
comes more sparse. It is interesting to note that a deeply 
bound state (e « —15 MeV) appears for I = (not shown 
in Fig. nj. This state corresponds to the Is level of the 
8 B nucleus. Since this state is occupied it is removed 
from our model space. We note also that some eigen- 
values may lie at very high energies. In principle, one 
could include these states in the system of the coupled 
equations. In practice, however, states at very high en- 
ergies are very weakly coupled to the ground state and, 
hence, do not influence the dynamics of the reaction. Be- 
sides, the inclusion of these states makes the calculation 
computationally more intensive and may even lead to 
convergence problems. For these reasons, in our calcula- 
tions those states above a certain excitation energy are 
completely removed from the coupled equations. This 
maximum energy will depend, indeed, on the particular 
reaction. In this case, we found that states above 10 
MeV have no effect whatsoever on the reaction observ- 
ables we are studying and, therefore, only those states 
with excitation energies below this value are retained in 
our calculations. Notice that states above 10 MeV have 
been also omitted from Fig. ^ 

In Fig. we plot the modulus of some components of 
the breakup S— matrix for a total angular momentum 
,/ = 150 which, assuming a classical trajectory, corre- 
sponds to a scattering angle of 10°. As noted in [20] at 
these angles breakup is mostly due to the Coulomb inter- 
action. In the standard CDCC, the breakup S— matrix 
is obtained by dividing the discrete S— matrix to the 
continuum bins by the square root of the bin width, 
Aki. In the CDCC-THO method, one could apply a 
similar procedure by assigning a width to each pseu- 
dostate. In 0] we used this approach to calculate the 
differential breakup cross section from the cross section 
to individual pseudostates. In particular, we assumed 
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FIG. 2: Breakup S-matrix elements for the total angular 
momentum J = 150 as a function of the asymptotic 7 Be+p 
relative momentum. 



that the width of the ith pseudostate is approximately 
given by Aj = (ei+i — ei_i)/2. In this work, we adopt 
a more sophisticated approach, previously proposed in 
in the context of the Gaussian expansion pseudostate 
method. In this method, the breakup S— matrix elements 
S a ':a(k), which depend on the continuous variable fc, as 
well as on the initial and final angular momenta, are ob- 
tained by an appropriate superposition of the discrete 
S— matrix elements S a >- a (ki) resultin g fro m the solution 
of the coupled-channels equations, as j2J| 
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S a '- a (k) ~ ^ ^ { 



:a{ki) 



(4) 



where (j^\(r) is the true scattering wave function for 

the partial wave £ and energy e = h 2 k 2 /2/i. The sum 
runs over the set pseudostates included in the coupled- 
channels calculation. The indexes a and a' denote the 
initial and final channels, that is, a = {gs; Lo,£o, J} 
and a' = {i;L,£, J} , where Lq (L) is the initial (final) 
projectile-target orbital angular momentum. 

The histogram represented in Fig. corresponds to 
the standard CDCC calculations, and the lines to the 
CDCC-THO results calculated with a THO basis with 
N P =N S =20 states. We verified that increasing the 
number of basis states does not change the calculated 
S— matrices, thus indicating the convergence of the THO 
method with respect to the basis size. Furthermore, the 
THO calculation with a N p =N s = 10 states gives already 
a result very close to result presented in Fig. |2j The fast 
convergence of the THO in this case can be attributed 
to the fact that the Coulomb interaction tends to pop- 
ulate mainly states in the continuum at low excitation 
energies, where the density of THO states is higher. 

In Fig. El we represent the elastic (upper panel) and 
breakup (lower panel) angular distributions. The stan- 
dard CDCC calculation is represented in both panels by 
a thick solid line. The thick dashed lines represent the 
CDCC-THO calculation with N = 40 states which, as 
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FIG. 3: Elastic (upper) and breakup (bottom) angular distri- 
bution for the 8 B+ 58 Ni reaction at 25.8 MeV in a model space 
including s— and p— wave states for the 7 Be+p continuum. 
In both plots, the thick solid line represents the converged 
CDCC calculation while the thick dashed line corresponds to 
the CDCC-THO calculation. The same potential parameters 
and partial waves are used for both calculations. The dotted 
line in the upper panel is the pure cluster-folding calculation 
without continuum. The thin lines in the bottom panel are 
the calculations including only nuclear breakup. 



can be seen, is almost identical to the standard CDCC 
calculation. The dotted line in the upper panel corre- 
sponds to the elastic calculation omitting the coupling to 
the continuum. Since both the CDCC and CDCC-THO 
use the same ground state wave function and interac- 
tions, this calculation is identical in both methods. The 
thin lines in the bottom panel are the calculated breakup 
angular distributions including only nuclear breakup. 

Next, we study the breakup differential energy cross 
section, da/dE. In the standard CDCC method, a natu- 
ral way to approximate this quantity is to divide the cross 
section for each final state by the bin width 0, El • In 
the THO method, we could proceed in a similar way, by 
defining a width for each pseudostate. Alternatively, one 
can calculate the breakup cross section from the continu- 
ous breakup S— matrices evaluated by means of Eq. iQJ. 
Both methods are compared below. 

In Fig.Q] we depict the differential energy cross section 
to £ — 1 continuum states obtained within the CDCC- 
THO method with several basis sizes, as indicated by 
the labels. The lines were obtained using the continu- 
ous S'-matrices given by Eq. <@J for different basis di- 
mensions. The circles represent the calculation with 
N s = N p = 60 states in which da/dE is approximated by 
dividing the cross section of the ith pseudostate by the 
width Aj = (ej+i — ej_x)/2. We see that both procedures 
to obtain da/dE give similar results. From this figure we 
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FIG. 4: Breakup cross section to p states as a function of the 
8 B internal excitation energy and integrated from zero to 180 
degrees. The lines correspond to the CDCC-THO calculations 
with different basis dimensions, evaluated from continuous S- 
matrices according to Eq. 0J. The solid circles represent the 
THO calculation with N s = N p = 60 states, in which the 
cross section to each pseudostate is divided by the estimated 
energy width (see text). 



can also infer that the calculation converges as the num- 
ber of basis states is increased. We notice that, at low 
excitation energies, the method converges even with a 
small basis. From this result we can conclude that the 
method is very suitable to describe Coulomb breakup re- 
actions. Convergence of the breakup cross section for ex- 
citation energies above 1 MeV requires a relatively large 
THO basis. This slower convergence rate is due to the 
fact that the density of THO states diminishes as the ex- 
citation energy increases, and more states are required to 
describe the excitation into this region. As the number of 
basis states increases, the solution of the coupled equa- 
tions becomes slower. However, we note that the actual 
number of states included in the coupled equations is re- 
duced with respect to the initial basis dimension. As ex- 
plained above, after the diagonalization only those states 
with excitation energy below 10 MeV are included in the 
coupled-channels calculation. In addition, the number of 
basis states is further reduced by removing those states 
with excitation energy below ~ 0.1 MeV. Once these two 
restrictions are applied, the number of continuum states 
actually included in the coupled-channels calculations for 
the case N = 60 is A^oup = 27 for both the s and p 
waves. We studied also the convergence of the method 
for the breakup to s-states. We found that, in this case, 
convergence is achieved with a smaller number of states, 
compared to the p-wave breakup. 

Finally, in Fig. we compare the CDCC (squares) and 
CDCC-THO (solid lines) breakup energy distributions. 
The latter corresponds to a basis with N s = N p = 60 
states for which, as we have shown above, good con- 
vergence of the THO method is achieved. We present 
separately the contribution of the I = (upper panel) 
and I = 1 (bottom panel) continuum states. We see that 
both energy distributions are in good agreement with the 
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FIG. 5: Breakup cross section as a function of the 8 B in- 
ternal excitation energy for the 8 B+ 58 Ni reaction at 25.8 
MeV integrated from zero to 180 degrees. The squares corre- 
spond to the standard CDCC calculations, whereas the solid 
lines represent the CDCC-THO calculation with a basis with 
N a = N p = 60 states. The separated contribution of s (upper 
panel) and p (lower panel) waves is displayed. 



CDCC calculation. 

In conclusion, in this work we have presented an exten- 
sion of the THO method, formerly applied to describe the 
£ = deuteron continuum, to generate a discrete repre- 



sentation of the continuum of a two-body system for any 
partial wave. As in the original formulation, the only a 
priori prerequisite for the application of the method is 
the knowledge of the ground state wave function of the 
system, either analytically or numerically. A local scale 
transformation, s(r), is then defined such that it converts 
the ground state wave function of the system into the 
harmonic oscillator ground state. Then, the THO basis 
is obtained by multiplying the ground state wave func- 
tion by a set of Laguerre polynomials expressed in the 
variable s(r). This THO basis constitutes a discrete and 
finite representation of the continuum part of the spec- 
trum with the same inter-cluster angular momentum as 
the ground state. For other partial waves, the method is 
generalized in a straightforward way by simply applying 
the same scaling transformation. The method has been 
applied to the reaction 8 B+ 58 Ni at subcoulomb energies, 
showing a fast convergence with the respect to the num- 
ber of basis states. Our results agree very well with those 
obtained with the standard CDCC method. In a recent 
work, we have successfully extended the method to de- 
scribe the continuum of a three-body system and the 
application of the method to describe the collision of a 
three-body loosely bound system by a target is underway. 
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